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Abstract 

Finite dimensional representations of extended Weyl-Heisenberg algebra are studied both from 
mathematical and applied viewpoints. They are used to define unitary phase operator and the 
corresponding eigenstates (phase states). It is also shown that the unitary depolarizers can be 
constructed in a general setting in terms of phase operators. Generation of generalized Bell states 
using the phase operator is presented and their expressions in terms of the elements of mutually 
unbiased bases are given. 
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1 Introduction 



The existence of phase operator, for quantized single mode electromagnetic field, was first postulated 
by Dirac pQ. However, since the harmonic oscillator spectrum is infinite, the creation and annihila- 
tion operators do not admit a naive polar decomposition. The infinite dimensional character of the 
harmonic oscillator Fock space constitutes a drawback in defining a phase operator in a consistent 
wa y OH]. To overcome this problem, Pegg and Barnett considered a truncated oscillator of finite 
dimension and defined the Hermitian phase operator [5] . But the truncated oscillator algebra has also 
one disadvantage from the algebraic point of view: its operators do not form a closed algebra. More- 
over, the truncated oscillator is not the only way to realize the Pegg-Barnett phase operator. Indeed, 
many approaches were proposed in the literature as for instance the refined version of phase operator 
proposed in [6J. This approach is based on the generalized Weyl-Heisenberg algebra which provides 
finite dimensional representations and where the hermiticity of the phase operator is automatically 
ensured. 

In Section 1, we consider a particular, but nevertheless very large class of extended Weyl-Heisenberg 
algebra. It comprises many well-known oscillator algebras, as particular or limiting cases, constructed 
in [H El El QUI EI]- We shall pay attention to extended Weyl-Heisenberg algebras possessing finite 
dimensional representations. Based on polar decomposition of raising and lowering generators, we 
construct algebraically the Hermitian phase operator. We discuss the temporal stability of the phase 
states (the eigenstates of the phase operator) and particularly the key role of the phase parameter 
ensuring this temporal stability. We study some quantum properties of the phase operator and compare 
with the results obtained using the Pegg-Barnett formalism. 

Section 2 of this work deals with the connection between the unitary phase operators and the 
so-called unitary depolarizer operators. This is mainly motivated by the fact that for finite dimen- 
sional systems, unitary operations have a prominent importance in quantum information processing. 
In particular, unitary depolarizer operators are useful for quantum information processing as for in- 
stance quantum teleportation, quantum dense coding |12[ [T3l Q3] as well as in investigating quantum 
optical systems. Note also that the notion of unitary depolarizer operators is deeply related to the 
entanglement phenomenon which has attracted a considerable attention as a crucial resource for quan- 
tum information processing involving higher dimensional systems. Here, we show clearly that unitary 
depolarizer operators can be constructed from the phase operator using the algebraic structures of the 
extended Weyl-Heisenberg algebra. We stress that our construction generalizes one presented recently 
in [15]. 

In Section 3, we show that the generalized Bell states introduced in |16] can be generated by means 
the generalized phase operator. It is important to mention that many efforts have been devoted, during 
the last years, to such states to understand bipartite entanglement in higher dimensional systems 
|17[ \18\ [T9l [20l 121) . The generalized Bell states provide a basis of maximally entangled states in 
a finite dimensional Hilbert space. In this respect, we examine interesting entanglement properties 
of some particular superpositions of generalized Bell states. Finally, since mutually unbiased bases 
[221 [23] can be derived from the phase states [6] and the generalized Bell states are defined by means 
of phase operator, it is natural to ask about the connection between generalized Bell states and 
vectors of mutually unbiased bases. In this sense, inspired by the results recently presented in [24] , we 
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derive such a relation taking advantage of the relation, between phase states and mutually unbiased 
bases, derived in [BJ. This establishes a useful link between mutually unbiased vectors containing 
single particle information and Bell vectors reflecting the maximal bipartite entanglement. Concluding 
remarks close this Letter. 

2 Extended Weyl-Heisenbeg algebra and phase states 

In this section we first introduce the extended Weyl-Heisenberg algebra. We indicate the constraints 
of the structure relations to get finite dimensional representations suitable to ensure the unitarity of 
the phase operator. Then, we construct the temporally stable phase states. The expectation value of 
the Hermitian phase operator is explicitly derived. 

2.1 Finite Fock space for extended Weyl-Heisenberg algebra 

The extended Weyl-Heisenberg algebra is an associative algebra over the complex field C, with gener- 
ators {a + ,a~,N} and the unit I , satisfying the structure relations 

[N,a~]=-a-, [N,a + ] = +a + , [aT ,a + ] = G(N), (1) 

where the non- negative function G(N) = (G(N))^ is some Hermitian analytic function of the number 
operator N. It is clear that for G(N) = I, one recovers the usual harmonic oscillator algebra. Also, 
this algebra reproduces the Wk algebra discussed in [TT] in the context of fractional supersymmetric 
quantum mechanics of order k and covers the extended harmonic oscillator worked out in [8, 9j. 

Let us consider the abstract Fock representation of this algebra by means of a complete set of 
orthonormal states {\n),n E N} which are eigenstates of the number operator N, N\n) = n\n). In 
this representation, the vacuum state is defined as a~|0) = and the orthonormalized eigenstates are 
constructed by successive applications of the creation operator a + . Indeed, we define the actions of 
creation and annihilation operators as 

a~\n) = y^n)e i[F{n) - F{n - 1)]v \n - 1), a+\n) = ^ F(n + l) e - t[F(n+1) - F{n) ^\n + 1) (2) 

where the structure function F(.) is an analytic function and should satisfy the conditions 

F(0) = and F(n) > 0, n = 1, . . . . (3) 

In what follows we shall denote the Fock space as T . On this space the operators a + and a~ are 
mutually adjoint, a + = (a - )^. It is easy to check that the structure function satisfies the 
following recursion formula 

F(n + l)-F(n) = G(n), (4) 

which gives by simple iteration 

n-l 

F(n) = Y, G(m). (5) 

m=0 
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We restrict ourselves here to generalized oscillator algebra defined through structure functions obeying 
the condition 



F(2s + 1) = 0, (6) 

where 2s is a positive integer. We assume that F{n) does not have zeros at non-negative integers values 
of n, except and 2s + 1. It follows that, in this case, the creation-annihilation operators satisfy the 
nilpotency relations (a~) 2s+1 = (a + ) 2s+1 = 0. This means that the corresponding representation is 
(2s + l)-dimensional. It is interesting to note that using ([5]), the condition © can be written as 

Tr G = (7) 

where the trace is over the (2s + l)-dimensional Fock space. The structure function F(N) can be 
factorized as 

F(N) = y s (2s + l-N) f(N) (8) 

where f(N) does not have vanishing eigenvalues for < n < 2s. Obviously, this algebra covers the 
usual harmonic oscillator for f(N) = I when s is large. It also includes the finite dimensional oscillator 
algebra A K (k < 0) defined in [6] when f(N) = I and k = — l/2s. Of course, as evoked above, the 
algebra introduced here covers many other variants of generalized harmonic oscillators, but we quote 
only the ones possessing finite dimensional representations which are interesting in the context of 
quantum optics. 

The algebra ([1]) defined by means of the structure function F(N) of the form ([8]) can be realized in 
terms of the usual single photon operator. For that end, one can write the creation and annihilation 
operators as 

a ~ = b ~ v^-^^W a+ = v/ /M\/ 1 -^r b+ ( 9 ) 

where b~ and b + and the number operator iV = b + b~ are the usual harmonic oscillator operators. As 
example of the function f(N), one may choose the function 

f(N) = (b-) m (b + ) m meN 

introduced in [7] to construct the multiphoton realization of generalized oscillator algebra. For our 
purpose the explicit form of the function f(N) is unessential through this work (except in the last 
subsection of the section 3 where we set f(N) = I). 

Using the algebraic structure of the generalized oscillator algebra, one can introduce an operator 
which generalizes the Hamiltonian a + a~ for the one-dimensional harmonic oscillator. Starting from 



a + a -\n) = F(n)\n) => H{N) = F(N) = a + a~. (10) 

We refer H(N) to as an Hamiltonian associated with the extended Weyl-Heisenberg algebra. The 
eigenvalue equation 

Tl 

H(N)\n) = F(n)\n) = — (2s + 1 - n) f(n)\n) (11) 
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gives the energies of a quantum dynamical system described by the Hamiltonian operator H(N). It 
is obvious that for f(N) = I and s — > oo, the Hamiltonian reads 

oo 

H(N)=Y^n\n)(n\, (12) 

n=0 

and for f(N) = I and s finite, one has 

2s 

H ( N ) = H £(2* + l-n)|n>R (13) 

n=0 

which coincides with the Hamiltonian introduced in [6]. In this particular case, it is interesting to 
stress that the ladder operators a + and a~ can be related to Stokes operators introduced in |25[ [26] 
to define a unitary operator representing the exponential of the phase difference between two modes 
of the electromagnetic field. This relation may be expressed as 

s+ = y/s a + S_ = y/s a_ s 3 = ^(N - s) 

and one can simply verify that the Stokes generators s+, s_ S3 close the su(2) algebra. It follows that, 
in the special case where f(N) = N(2s + 1 — N)/s, the extended Weyl-Heisenberg generators can be 
realized a la Holstein-Primakoff by means of two independent ordinary bosons (two electromagnetic 
field modes). Remark that the parameter s can be viewed as controlling the deformation scheme. 
Indeed, when s becomes large one has the standard harmonic oscillator algebra. We mention that the 
idea of extended Weyl-Heisenberg is mainly inspired by the polynomial deformation of Lie algebras 
introduced in [271 I2H] and extensively discussed in the context of quantum algebras [29j [301 EI]- 
Finally, we emphasize that quantum systems described by a nonlinear Hamiltonian H(N) are familiar 
in the context of nonlinear quantum optics as for instance electromagnetic field propagating trough a 
nonlinear Kerr medium [32J. 

2.2 Temporally stable phase states 

The Hilbert space T is (2s + l)-dimensional. The actions of a~ and a + on T are given by the equation 
([2]) supplemented by 

a + |2s)=0, (14) 

which easily follows from the calculation of (2s|a _ a + |2s). 

Let us consider the polar decomposition of the creation a + and annihilation a~ operators: 

a~ = E\/F(N) ^ a + = y/F(N) (£) f . (15) 

It follows that the operator E satisfies 

E\n) = e'FW-^-^ln - 1) (16) 

for n = 1, 2, . . . , 2s. For n = 0, the action of E is defined by 

E\Q) = e i\m-FWA>p\2a) (17) 
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so that (|16|) is valid modulo 2s + 1. It follows that we have 

(£) f \n) = e-^ F ^- F ^f\n + 1), (18) 

where n + 1 should be understood modulo 2s + 1. It is important to stress that the operator E 
is unitary. Therefore, equation (|15p constitutes a polar decomposition of creation and annihilation 
operators a + and a~. 

To construct the phase states associated with the finite dimensional algebra under consideration, 
let us look to the eigenstates of the operator E. For this, one has to solve the eigenvalue equation 

2s 

E\z)=z\z), \z) = ^C n z n \n) (19) 

n=0 

with z G C. According to the method developed in [6], it is easy to see that the eigenvalues z should 
satisfy the discretization condition 

z 2s+1 = 1, (20) 
and subsequently the complex variable z is a root of unity given by 

z = q m m = 0,1,... ,2s, (21) 

where 

g:=e *rf/(2«+i). (22) 

Then, the normalized eigenstates \z) = \m, ip) of E are given by 

1 2s 

\m, if) = - V e- iF{n ^q mn \n). (23) 

v n=0 

The states \m, (/?), labeled by the parameters m G Z/(2s + 1)Z and i^GK, satisfy 

27T 

S|m,¥>) = e jem |m,^) ^ m = m— — , (24) 

As + 1 

which reflects that i£ is indeed a phase operator. In the particular case (p = 0, the states \m, 0) 
correspond to an ordinary discrete Fourier transform of the basis {|n) : n = 0,1,..., 2s} of the 
(2s + l)-dimensional space T . At this level, we would like to emphasize the key role of the parameter 
ip introduced from the beginning in the definition of actions of creation and annihilation operators 
([2]). First, it ensures the temporal stability of the phase states \m,ip) under "time evolution": 



e 



- iH ^ t \m,ip) = \m,(p + t) (25) 



for any value of the real parameter t. Secondly, if one ignores the parameter <p,i.e. <p = 0, the 
phase states ([23]) reduce to ones derived by Pegg and Barnett using the truncated harmonic oscillator 
[5]. This means that for two different extended Weyl-Heisenberg algebras characterized by different 
structure functions, we will end up with the same phase states. Then to differentiate between them, 
the phase parameter is necessary. In the same spirit, it must be noticed that the 577(2) phase states 
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obtained in [33J are similar to ones obtained by Pegg and Barnett despite the fact that the involved 
symmetries are different. Hence, to avoid such a problem, we introduced the parameter (p which allows 
us to keep the trace of the symmetry and the dynamics of the system under consideration. Note also 
that the parameter p plays a key role in relating phase states and mutually unbiased bases [6]. The 
phase states have the following remarkable properties. For fixed (p, they satisfy the equiprobability 
relation 

\{n\m,ip)\ = - 1 n,raeZ/(2s + l)Z. (26) 

y2s + 1 

They constitute an orthonormal set, for fixed ip, 

(m, p\m' , ip) = 5 m ^ m > m,m' £ Z/(2s + 1)Z (27) 
and satisfy the closure property 

2s 

^2\m,<p)(m,<p\=I. (28) 

m=0 



Finally, the overlapping between two phase states \m',(p') and \m,(p) is given by 

2s 

2 s + ~ 



71=0 

where 



1 2s 

<m^|mV> = 27TT^ gP(m ~ m ,iP ~ V,n) ' (29) 



2s + 1 

p(m — m! , ip — if' , n) = — (m — m')n H (ip — ip')F(n) (30) 

27T 



and q is defined in (|22p . This shows that the temporally stable phase states are not all orthogonal. 
2.3 Expectation value of the phase operator 

In the previous subsection, we defined the unitary phase operator and the phase states. We shall 
calculate the expectation value of the Hermitian phase operator defined by 

E = e ie . (31) 

It is clear that the Hermitian phase operator can be written as 

2s 

6 = m \m, ip) (m, ip\ (32) 

m=0 

where the summation is modulo 2s + 1. Using the expression of the phase states (|23p . one obtains the 
expansion of the Hermitian phase operator, in the Fock states basis, as 

= 7^2 E E e-^-W], eM | n)(n /| (33 ) 

^ ^ > n=0n'=0 

where the objects 6(n,n') are given by 

2s+l 



9(n,n') = E m q m{n ~ n>) = (2s + 1) 



ra=l 



gU—n 

S $n,n' + z Z—7 (1 - 8 V 



(34) 
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Using this result, we are now in position to compute the phase-number commutator. In the number 
basis, it can be expressed as 

2s 



[e,N] 



2tt 



2s + 



in — n 



1 



-i[F(n)-F(n')}<p 



\n){n'\ 



(35) 



This equation shows that the diagonal matrix elements of the phase-number commutator vanish and 
coincide with the result obtained by Pegg and Barnett [5] when <p = 0. This indicates clearly that the 
generalized Weyl-Heisenberg algebra gives an alternative way to deal with finite dimensional harmonic 
oscillator in investigating the phase properties of an electromagnetic field. In this respect, the form of 
the phase states suggests us to consider a normalized pure state of the form 

2s 



II } 



(36) 



n=0 



to evaluate the phase probability distribution and the phase expectation value. The parameter a is 
real and the coefficients $ n are real and positive satisfying 

2s 
n=0 

The vector |$) is similar to the so-called partial phase state introduced in [5) (up to phase factor 
involving <p) and the phase probability distribution is 



{m,ip\$)\ 2 = ——- + —— $n®n>cos[(n-n r )(a-e m )}. 



n>n' 



It is bounded as follows 
1 



2s + 1 



2s 



n=0 



Is 
^ n=0 



(37) 



(38) 



Using ([33]) . the expectation value of the Hermitian phase operator reads 

„ 2s 2s 

2tt 



(6) = ($|0|$) = (2g + 1)2 EE ®n$n>e ia{n - n,) (n,n') 



n=0n'=0 



which can be written as a sum of two terms 

(6) = (G) d + (G) nd 



(39) 



(40) 



where (@)d stands for diagonal contributions (n = n') and (0) n d is the contribution from the off- 
diagonal terms. These quantities are given by 

2tts 



(©}c 



and 



(G) 



2tt 



nd 



2s + 



I E 



2s + l' 
sin((a- 2jri ){n-n')) 



Sin (2^l( n - n ')) 



(41) 



(42) 



n>n' - \2s-fl 

It is clear that for a = n/(2s + 1), the off-diagonal contribution vanishes and the Hermitian phase 
operator expectation value is 2ns /{2s + 1). 
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3 Unitary depolarizer operators and generalized Bell states 



As mentioned in the introduction the extended Weyl-Heisenberg algebra not only introduces math- 
ematical tools for the analysis of finite quantum systems but also can be applied, for instance, to the 
context in quantum information processing via the so-called the unitary depolarizer operators and the 
generalized Bell states. In this section, we shall present the derivation of unitary depolarizer by means 
of phase operator. We also obtain generalized Bell states via the action of unitary phase operator and 
establish a relation between them and vectors generating mutually unbiased bases (MUB). 

3.1 Unitary depolarizers operators 

In a d-dimensional Hilbert space, the unitary depolarizers on a domain T> are the elements of set 

S d = {U a | UaPl = UlU a = I, a e V}, (43) 

which satisfy the relation 

\ U a OUt = (TrO) I (44) 

for any operator O defined on the Hilbert space. The phase operator associated with extended Weyl- 
Heisenberg algebras discussed in the previous section offers a simple way to define unitary depolarizers. 
For this, we rewrite the unitary phase operator E as 

2s 

E = Y^ e i{ - Fi > n+1 ^ F ^\n){n + 1|, (45) 

n=0 

from which one can obtain 

2s 

E k = e i{F{n+k) - F( - n))tp \n)(n + k\ modulo (2s + 1). (46) 

n=0 

The phase states ([23]) are eigenstates of the operator E k . Indeed, one has 

E k \m,ip) = e ik9m \m,p). (47) 

Any operator O, acting on the Fock space T, can be expanded in the phase states basis {\m,(p),m G 
Z/(2s + 1)Z} as follows 

2s 2s 

0=^2 Yl °m,mf\m,(p)(m',(p\. (48) 

m=0 m'=0 

It is simple to check that 

2s 2s 

Y^E k OE kj = (2s + l)J2 m,m\m, V ){m, l p\. (49) 

k=0 m=0 

Next, we define the unitary operator 

V(<p) = e^W+^f (50) 
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in terms of the number operator and the structure function F(.). Using the equations (|49p - (|5L)p . the 
expression of the phase states (|23p and performing integration over (p, one gets 



- [ +W dipJ2v{(p)E k OE k W(<py = TiO I. (51) 



This constitutes a generalization of the definition of unitary depolarizers operators discussed in [15j . 
Then, the definition of unitary depolarizers (eqs. (|4"3j) and (|44p) can be extended to unitary operators 
labeled by continuous parameters belonging to some compact domain. It is important to note that 
for quantized values of the parameter ip, given by 

2ir 

<p = ——l, i€Z/(2« + l)Z, 
2s + 1 



one can verify that the discrete analogue of the equation (I5ip is 

2s 2s 

^2 ^2 ViE h OE k ^Vi = {2s + l)TrO I (52) 

k=0 1=0 

where Vi = V(2ttI/(2s + 1)). This agrees with the result obtained in |15| . So, the unitary depolarizer 
operators can be expressed in terms of the phase operator. To close this subsection, it is important 
to stress that the construction of unitary depolarizers is done without specifying the nature of the 
structure function F(.). This provides us with a general scheme to define the unitary depolarizers 
for a large class of extended Weyl-Heisenberg algebras and can be adapted to quantum systems with 
higher symmetries. 

3.2 Generalized Bell states and MUB 

In this last part of our work, we shall focus on extended Weyl-Heisenberg algebra defined by the 
structure function 

F(N) = ^-(2s + l-N). (53) 

In this case, we will show that the generalized Bell sates can be generated from the maximally entangled 
states using the phase operator. Recall that the generalized Bell states were first introduced in [16] 
to study quantum teleportation for higher dimensional quantum systems. We discuss briefly the 
entanglement properties of some particular superpositions of generalized Bell states and establish a 
connection between the generalized Bell states and the vectors of mutually unbiased bases. 
The maximally entangled state of two (2s + l)-dimensional systems is 

1 2s 

|MES) = -_=£|n)®|n>. (54) 

n=0 

The action of the operator E k I on the state |MES) leads to 

2s _ 2s 



E k 



I ; 1 V In) ® In) = ; 1 V e i(F(n+k)-F(n))<p\ n \ g, U + fc\ (55) 
n=0 n=0 
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Using (|53|) and setting ip = — |7r (this discretization condition of </3 is similar to one introduced 
in [6] to get mutually unbiased states from phase states), one can write the equation (|55l) as 

1 23 

E k I -=_ £ |n) ® |n> = g" (fc - 2s - 1} |^, p ), (56) 

n=0 



where the states \^?k,p) are the generalized Bell states introduced in [16]. They are given by 

1 



■,P 

1 23 

qr k ) = ^=J2q n P\n)®\n + k), (57) 

It follows that the generalized Bell states can be obtained from the maximally entangled states (|54p 
by means of a successive application of the unitary phase operator. 

Finally, it is interesting to note some characteristics of some particular superpositions of generalized 
Bell states ») an d their expressions in terms of vectors belonging to mutually unbiased bases. 
First, it is remarkable that the superposition of the (2s + l) 2 maximally entangled Bell states 

k=0 p=0 v n=0 

is a completely separable state. In other words, the superposed state is disentangled despite that 
it is a superposition of maximally entangled states. This is not always true. Indeed, the following 
superposition 

2s 2s 2s 

£ = -== 5>>^£ q nk \n + k). (59) 

k=0 V ZS + 1 n=0 k=0 

gives a maximally entangled state. The question concerning the entanglement of superposed Bell 
states is very rich and deserves more investigation. We hope to report on this issue in a forthcoming 
work. There is also a nice connection between the superposition of Bell states and the vectors of 
mutually unbiased bases B p introduced in [6] as 

1 2s 

B p = {|C> = -7^=j E q* n{2s+1 - n) q nm \n), m G Z/(2« + 1)1} p G Z/(2a + 1)Z. (60) 

n=0 

Indeed, the superposed states (|58|) can be expressed in terms of the mutually unbiased vectors \<$n). 
This expression is 

2s 2s 

L 

2s 



k=0 p=0 

and the state (|59p can be written as 

2s 2s 



rViEEi^) = i )®i^) 5 (si) 



= X^-" |n>®|^>. (62) 

fc=0 n=0 
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The vectors 



involved in the above superposition belong to the basis Bq. One can go further to 



find superpositions of maximally entangled states related to vectors with p ^ 0. Indeed, it is easy 
to see that the action of the unitary operator I ® e lF ( N ) s ^ with 5 — 1 " 



fJV(2s+l-iV) 



2s 
fc=0 



k.k) 



2.s 



n=0 



n) (g> 



, , ^jwr, on the states (|62j) gives 

(63) 



This establishes a simple relation between the generalized Bell states and the vectors of the mutually 
unbiased bases. To close this section, we notice that the vectors of mutually unbiased bases B p can 
be defined as the quadratic Fourier transform of the Hilbert space basis [Mj. Remark also that the 
generators of generalized Pauli group, which constitutes an important ingredient to deal with the 
relation between generalized Bell states and MUBs [24], can be defined as in 



4 Concluding remarks 

This Letter is based on the extended Weyl-Heisenberg algebra. This algebra is interesting in many 
respects. It describes in a unified way a large class of generalized harmonic oscillators encountered 
in the literature. The extended Weyl-Heisenberg algebra, under some assumptions, possess finite di- 
mensional representations. The nonlinear structure function characterizing the algebra can take into 
account some nonlinear effects that may occur in the quantum description of quantized modes of the 
electromagnetic field [35]. In connection with the extended Weyl-Heisenberg algebra, this work deals 
with the construction of a unitary phase operator and the determination of its temporally stable eigen- 
states (the so-called phase states) . This construction is valid for a large class of generalized harmonic 
oscillators. We also discussed the derivation of unitary depolarizer operators from the unitary phase 
operators. We have shown that the formalism, presented in [15], using the Pegg-Barnett operator can 
be generalized to a large variety of extended harmonic oscillator. Finally, we discussed the generation 
of generalized Bell states [16] using the phase operator. We also investigated the degree of entangle- 
ment of superpositions of Bell states and established a connection between them and the vectors of 
mutually unbiased bases. 
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